Singlet and triplet BCS pairs in a gas of two-species fermionic polar molecules 
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We investigate the BCS pairing in a mixture of fermionic polar molecules with two different 
hyperfine states. We derive a set of coupled gap equations and find that this system supports both 
spin-singlet and -triplet BCS pairs. We also calculate the critical temperatures and the angular 
dependence of order parameters. In addition, by tuning short-range interaction between inter- 
species molecules, the transition between singlet and triplet paired states may be realized. 
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Introduction. — There have been growing interests 
in studying the properties of dipolar Fermi gases since 
the experimental success in making high phase-space- 
density fermionic KRb gas [U, Of particular inter- 
est is the anisotropy of electric dipole-dipole interaction 
(EDDI), whose partially attractive character provides a 
mechanism to form anisotropic BCS pair. You and Mari- 
nescu 3 first realized that the p-wave paired BCS states 
could be achieved for fermionic atoms inside an external 
dc field. Baranov et al. further showed that the order 
parameter of a BCS pair is the superposition of all odd 
partial waves Subsequently, the critical temperature 
of the superfluid transition and its relation to the trap ge- 
ometry were investigated [4]. In those studies, the Fermi 
surface (FS) of system was assumed to be a sphere. How- 
ever, it was found via variational method that the Fock 
exchange interaction deformed the FS a| . This result was 
also confirmed by both numerical 0, [? | and perturbation 
calculations d, Q . Very recently, the role played by Fock 
exchange interaction in BCS pairing is studied 0. 

All above mentioned theoretical work focus on the 
spinless dipolar fermions. Adding the spin degree of 
freedom, ultracold gas of polar molecules may offer even 
richer physics. When loaded into an optical lattice, spin- 
i polar molecules can be used as a toolbox for lattice 
spin models p^ . Furthermore, it was shown that a two- 
component dipolar Fermi gas may have a ferro-nematic 
ground state [11[. 

In the present work, we study the BCS pairing in a 
mixture of fermionic polar molecules with two differ- 
ent hyperfine states, cr =t and [. The electric dipole 
moments d of all molecules are orientated along z-axis 
such that the spin-independent EDDI becomes Vd(r) = 

CdT"'^ - with Cd = 4^^. In addition, we as- 
sume that inter-species molecules also interact via short- 
range interaction Vb(r) = Co5{v), where Co — 
with Us being the s-wave scattering length and m the 
mass of the molecule. The total interaction potential 
V = Vd + Vq conserves the number of particles in indi- 
vidual spin component. We show that novel spin-singlet 
even partial wave and spin-triplet odd partial wave pairs 
can appear in this system. 

Gap equation. — We consider a homogeneous gas of 



two-species fermionic polar molecules with total number 
density 2n. For simplicity, we assume that n^- = n ^ = n. 
In momentum space, the second quantized Hamiltonian 
reads 
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where fj, — — fii is the chemical potential, V is the 
volume of the system, and V = Vd -\- Vq is the Fourier 
transform of the interaction potential. One should keep 
in mind that, in Eq. ([1]), contact interaction only exists 
between t and | molecules. 

Under Hartree-Fork approximation, Eq. ^ becomes 
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where the dispersion relation is eko- = — sp — ^a (k) 

with sf = fJ- — nCa = ^{Qir'^nY/^ being the Fermi 

energy and Scr(k) = V"^ Z]q^d(k - q)/qcr the self en- 
ergy. To the lowest order, we choose the Fermi occu- 
pation number /qo- = (c^o-Cqo-) to be that at zero tem- 
perature. As a result, dispersion relations become spin 
independent, denoted as ek- The Fock exchange interac- 
tion deforms FS such that, in the vicinity of FS, Ek can 
be approximately expressed as Q 
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where a is the deformation parameter, which can be 
found using either variational or numerical method [B, 
Q. Moreover, at the weak interaction limit 

D = nCd/ep « 1, (4) 

the perturbation calculation yields a = 1— @, [1], 
where D dimensionless dipole interaction strength. We 
point out that, for all results presented in this work, a is 
obtained variationally. 
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The order parameters for the BCS states in the Hamil- 
tonian ^ are defined as 



(5) 



Using either Bogohubov transformation J^] or the stan- 
dard Green's function method [l^, we obtain 



2^ tanhl ) , (6) 



where £'k<j = ye^ + I^CTc(k)p is the energy of single- 
particle excitation and /3 = l/(fcBT) is the inverse tem- 
perature. Substituting Eq. ^ into ([5]), one find a set of 
self-consistent equations for order parameters 
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We note that the order parameters corresponding to dif- 
ferent BCS paired states are coupled in above equations. 
Similar gap equations also appear in superfiuid ^He sys- 
tem [l^. The integral in Eq. ([7]) formally diverges. To 
regularize the interaction, we follow the standard renor- 
malization procedure 0], which, to the first Born approx- 
imation, leads to the gap equations 
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We emphasize that these equations are only valid in weak 
interaction regime Eq. as we have ignored the higher 
order contributions from the interaction. 

Symmetries of the order parameters. — Up to now, we 

Att A|x 



have expressed the order parameters A = 



AiT Aa 



in the uncoupled spin space. While in the basis of coupled 
spin {\SAI)}, one will have singlet and triplet pairs corre- 
sponding to total spin 5* = and 1, respectively. Singlet 
pairing requires the gap parameter to be antisymmetric 

such that As — Asi(T2 — ( 'I ] , where cr2 is the 



, -As ^ 

second Pauli matrix. Obviously, As(— k) = As(k), and 
as we shall show, it is a superposition of even partial 

waves, i.e., As(k) = Ee.en i (k)Ylo(^- 

Following the convention of ^He 13], we may define 
a vector d(k) in spin space by combining three spin 
components {M = 0, ±1) of the triplet pair: At = 

(T^ are Pauli matrices. It can be easily seen that 
Atj = A|| = -di + id2, At-i = All = di + ic?2, 
and Af^o = A|| = A|| = ^3. In momentum space, 
At^M(k) are superpositions of odd partial waves, i.e., 

At!M(k) = Eodd;Ai^yfc)y,o(k). 



Spin-singlet pairing. 



For deformed FS, we intro- 



duce rescaled momentum k such that k^/y — a ^I'^k^jy 
and fc^ = ctkz. Moreover, we define = ?i^k^/(2M) — 
Ei?, which allows us to rewrite the order parameter as 
As(k) = As(^j^,njj) with — k/fc being an unit vec- 
tor. Using the fact that pairing is mainly contributed by 
states near FS, we introduce a characteristic energy cUs 
(of the order of ep) to single out the contribution from 
— ZUs < — ,4]. In weak coupling regime, we have 
ujs ^ |As|, ksTc'^K The value of ZU^ will be determined 
self-consistently. After some tedious calculations, we ob- 
tain from Eq. ^ that 

A.Kt.nt) = - J ^H'(T,n,)K.(f|„n|;;0,n5) 
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where W^I\n^)_^J^ d^^r]^Ha.nh[Tj^^/i2kBT)] and 
775. = Y^C^ + I As(0, nk)p. Furthermore, the integration 
kernel takes the form 
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Close to the critical temperature Tc , the gap goes 
to zero and rjq |^q|. For the first term in rhs of Eq. 
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©, we have C = W (t^^', rik] = In 

with 7 = 0.5772 being the Eulcr constant. In contrast 
the last two terms in Eq. ^ do not contain the large 

logarithm In (zJs/kBTc'^^^ , indicating that they are only 

important for the preexponential factor in the expression 

of critical temperature (J]. For the purpose of determin- 
es) 

ing Tc , they are neglected such that we find on FS 

As(0,nk) = -C J ^7^s(0,nk;0,nq)As(0,nq), (10) 

which is essentially an eigenvalue equation of integral op- 
erator with kernel 7?,s(0, rik; 0, riq). Finding the highest 
critical temperature is equivalent to finding the lowest 
negative eigenvalue of this eigenvalue equation. Even in 
case a 7^ 1, it can be shown that the eigenstates corre- 
sponding to lowest eigenvalue is independent of the az- 
imuthal angle of rik [4,]. This observation allows us to 
integrate over the azimuthal variable of Eq. (fTU|) to ob- 
tain an eigenvalue equation 



dyK.s{.x,y)(t){y) = X(t>{x), 



(11) 
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FIG. 1: The a dependence of As for kpan = (solid line) and 
At (dashed line). 



where x = cosflnj., y — cos6'n^, and lCs{x,y) = 
[Zla{x,y)\x - y\ - l] + kpas/in^D) with Zq, = 



-1/2 



and fcp = 



{Aa^ + (1 - a^)[(a; - y)^ - a^{x + j/)^]} 
\f2MeF jh being the Fermi wave vector. Clearly, the 
first term in /Cs originates from EDDI, while the second 
one corresponds to contact interaction. 

For singlet pairing, we focus on the subspace of the 
eigenstates with even parity, i.e., ^(a;) = (/>(— x). Assum- 
ing that \s is lowest eigenvalue of Eq. (fTTj) corresponding 
to an even eigenstate 0s (x) , the critical temperature can 
then be expressed as 



fc^T(^)^^exp 



is) 

In addition, at Tc , the order parameter on FS becomes 
As(0,nk) = Aso(?!)s (cos6'„j.), where 0s (cos 6'„^ ) repre- 
sents the angular dependence of the order parameter on 
FS. 

If the Fock exchange interaction is ignored, we have ex- 
actly <i)s{x) = M COS (^x^/3J\X^\j , where A/" is a normal- 
ization constant and |As| is the largest positive root of 

the equation t cos 1/37*+ (1 + ^ffs^ ) ^in ^/sjt = 0. 
It can be shown that As (< 0) is a decreasing function of 
kpOLs/D. For negative kpas, our result reduces to that 
for s-wave pairing when D = 0. In case a 7^ 1, As and 
(j>s can only be determined numerically. As shown in Fig. 
[1] As is an increasing of a for kpag = 0, indicating that, 
in the absence of contact interaction, the Fock exchange 
interaction enhances singlet BCS pairing. 

Next, we compute the characteristic energy Us ■ Follow- 
ing the Gor'kov and Melik-Barkhudarov's approach ^4, 
Il4l. we obtain 
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where the order parameter, at Tc , satisfies the equation 
As(&,nfc)^ 1 /'^^^(^_^„_.o,„_)^^(cos^q)(14) 
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FIG. 2: (color online) (a) The D dependence of Tc for 
kpas — (solid line) and —1 (dashed line), (b) ri"' as a 
function of kpas for D = 0.1 (solid line) and 0.15 (dashed 
line), (c) As (ftp, kz)/Aso at critical temperature for kpaa — 
and D = 0.01, where kp = {kl + fe^)^/^. 



ZUs can be obtained by substituting Eq. (|T4)) into (|13p . 
We remark that, in the rhs of Eq. the second term 

also contains a divergent term which cancels the diver- 
gence from the first term. In general, ZUs is a function of 
both D and kpas, indicating that it also implicitly de- 
pends on a. For simplicity, we use a = 1 to calculate Wg, 
which yields ZUs — 0.36eF for kpas — and ZUs — 0.54ei;- 
for D = 0. The latter is exactly the value for s-wave 
pairing. 

In Fig. [21 (a), we present the D dependence of the crit- 
ical temperature for singlet pairing. For kpUg = 0, only 
dipolar interaction is responsible for pairing, therefore, 

(s) . 

Tc is a monotonically increasing function of D. How- 
ever, for kpGs = —1, Tc^\d) becomes a concave curve. 
This result suggests that even though both attractive 
contact interaction and dipolar interaction contribute to 
pairing, they also compete with each other. In fact, the 
deformed FS makes it difficult for contact interaction to 
form s-wave pairs, which is responsible for the drop of 

(s) 

Tc at small D. When dipolar interaction dominates, 

(s) 

Tc becomes a increasing function of D again. We also 

(s) 

verify that if the FS deformation is ignored, Tc would 
always be increasing function of D. Figure [5] (b) shows 

( s) 

the kpGs dependence of Tc for given D. As a is fixed by 

(s) 

D, Tc always increases as one increases the strength of 
contact interaction. In addition, when short-range inter- 
action dominates, the one with smaller FS deformation 
(smaller D) will eventually has a higher critical tempera- 
ture. Finally, the angular dependence of order parameter 
for spin-singlet pair at critical temperature is presented 
in Fig.[2](c). 

Spin-triplet pairing. — At critical temperature, the 
order parameters for triplet pairs At,M are decoupled. 
In addition, they satisfy gap equations with exactly the 
same form. Therefore, we denote the triplet order pa- 
rameter as At. Since the critical temperature and order 
parameters can be obtained by following the same proce- 
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FIG. 3: (color online) (a) The D dependence of Ti*\ (b) 
At{kp,kz)/ Ato at critical temperature for D = 0.01 
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FIG. 4: The dipolar interaction strength dependence of the 
critical scattering length. 

dure as that for singlet pairing, here we shall only present 
our results. 

The critical temperature Ti'-* for triplet pairing 
can be expressed as /cbT]*^ — ^^^^ exp ^— ^^rjjyj^^ , 
where At is the lowest eigenvalue of the equation 
J^-^^dylCt{x,y)(l){y) — X(j){x) in the subspace of the 
eigenstates with odd parity. The integration kernel is 
lCt{x,y) = Ha{x,y) \x — y\ — 1, which is essentially K-s 
with kpcis = 0. Assuming that the corresponding eigen- 
state of A( is <j)t {x) , the triplet order parameter on FS is 
At(0,cos6'nj-) = AtQ(j)t{cos6ji^). The a dependence of At 
(Fig. [T]) indicates that the triplet pairing is also enhanced 
by the Fock exchange interaction. 

The characteristic energy ZUt and the triplet order pa- 
rameter At(fj(., njj.)/Ato at critical temperature satisfy 
similar equations as Eqs. and except for that 
TZs is now replaced by TZt which does not contain the 
contact interaction. In Fig. [3] (a), we plot the dipolar 
interaction strength dependence of the critical temper- 



ature Tc , where we have used ujt — 0A2eF, the value 
corresponding to a = 1 case. The typical order parame- 
ter for triplet pair is plotted in Fig. [3] (b) . We point out 
that in case the FS deformation is ignored, our results 
for triplet pairing reduce to those obtained by Baranov 
et al. [1]. 

Singlet versus triplet pairing. — For kpag — 0, it can 
be seen from Fig. [l]that we always |As| < |Af |; by further 
assuming a — 1, we also find that tJs < 'uJt- These re- 
sults suggest that, without attractive contact interaction, 
the critical temperature for singlet pairing is always lower 
than that for triplet pairing. However, for kpas < 0, ri'^ 
increases as one increases the strength of short-range in- 
teraction such that beyond a threshold kpal, we have 
ri'*^ > ri*^ The critical s-wave scattering length is de- 
termined by equation 7r_D In ""(^^"a) = _L _ — — 1 — 
where we have explicitly expressed uJs and A^ as func- 
tions of kptts- In Fig. [4j we present the D dependence 
of the critical scattering length. We point out that when 
ri"^ = Ti*\ the singlet and triplet pairs may coexist in 
the system. 

Conclusion. — We have studied the BCS pairing in 
a ultracold gas of two-species fermionic polar molecules. 
We show that there exist two types of BCS pairs: spin- 
singlet even partial wave and spin-triplet odd partial 
pairs. We calculate the critical temperatures and angu- 
lar dependence of the order parameters. It is found that, 
without attractive contact interaction, the critical tem- 
perature for triplet pairing is always higher than that 
for singlet pairing. The opposite could also be true if 
one tunes the scattering length to a large negative value. 
Compared to the spin-triplet p-wave pairing in superfluid 
•^He, the spin-triplet (-singlet) pairing studied here cou- 
ples all odd (even) partial waves. In addition, by carefully 
tuning the contact interaction, singlet and triplet pairs 
may coexist in this system. Armed with the rich ex- 
perimental control mechanisms developed in cold atom 
physics, two-component fermionic polar molecules may 
provide an ideal platform for exploring new physics be- 
yond what superfluid '^He could offer. 
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